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Abstract 

We investigate the distribution of instanton sizes in the framework of a simplified 
model for ensembles of instantons. This model takes into account the non-diluteness 
of instantons. The infrared problem for the integration over instanton sizes is dealt 
with in a self-consistent manner by approximating instanton interactions by a repul- 
sive hard core potential. This leads to a dynamical suppression of large instantons. 
The characteristic features of the instanton size distribution are studied by means of 
analytic and Monte Carlo methods. In one dimension exact results can be derived. 
In any dimension we find a power law behaviour for small sizes, consistent with the 
semi-classical results. At large instanton sizes the distribution decays exponentially. 
The results are compared with those from lattice simulations. 



1 Introduction 

1.1 Instantons in gauge theories 

In non-abelian gauge theories different topologically nontrivial configurations have been 
made responsible for non-perturbative features |], 0, [| f|, ||. An important class are 
instantons, which are solutions of the Euclidean field equations with non-vanishing topo- 
logical charge ||. They give contributions to the saddle-point approximation of Euclidean 
functional integrals, which lead to non-perturbative effects [[?], §]. For a review see |L0 . 



In the dilute gas approximation || one considers superpositions of single instantons as 
quasi saddle points of the action. These configurations are characterized by the instanton 
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positions {aj} in four-dimensional space-time, the instanton sizes {pj}, and other internal 
parameters. For a single instanton the action is 
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where g is the bare gauge coupling constant. Taking into account quadratic fluctuations 
around the one-instanton solution |7], || its contribution to the functional integral is 



Z 1 = Jd 4 & JdpCp- 5 exp^- 



57T 



's 2 (l/p)J' {) 

where g(p) is the running coupling. In the case of supersymmetric Yang-Mills theory this 
formula has even been established at the two-loop-level, in ordinary Yang-Mills theory 
there are higher-order corrections to the integrand ||. In one-loop order the running 
coupling obeys 

61og(/i/A) 

where 

■ t 

and A is the renormalization-group invariant scale parameter, so that 

Z 1 = Jd^a JdpCp~ 5 (pA) b . (5) 

Therefore the integrand is proportional to p b ~ 5 and increases with increasing instanton 
size. 

The space-time integral over the instanton position gives the usual volume factor, 
which is needed in the large volume limit in order to get an extensive free energy. The 
integral over p, however, represents an infrared problem. Where the instanton density 
becomes important, for p > 0(1/ A), we leave the region of valididty of the semiclassical 
expansion because the running coupling becomes too large. 

It should be noted that the apparent infrared divergence in (|5|) is an artifact of us- 
ing the one-loop formula for g 2 (l/p) long after it has become invalid, i.e. for p > 1/A. 
Nevertheless we are confronted with the infrared problem in the size integration for the 
single instanton contribution Z\. If the semiclassical approximation is meaningful at all, 
a solution of this problem in the context of the full instanton ensemble is required. 

The quasi saddle points composed of any number of instantons and anti-instantons are 
treated as independent in the dilute gas approximation. Consequently their contributions 
exponentiate in the usual way. The problem with the integration over the sizes persists 
and has to be dealt with. The simplest way is to cut the integrations off at some ad-hoc 
value p c . But since the integrand increases with increasing pj the dominant contribution 
comes from large pj near the cut-off where the assumption of diluteness fails. Moreover 
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the introduction of an ad-hoc cut-off leads to inconsistencies with the renormalization 
group 
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In order to solve the problem of the instanton size integrations is has been proposed 
that instanton sizes are cut off in a dynamical way [11, 12| . The dynamical cut-off 
should originate from configurations where instantons start to overlap. Configurations of 
overlapping instantons have an action which deviates from the sum of the single instanton 
actions. Therefore large instantons feel an interaction. Additionally, the fluctuations 
around the multi-instanton configurations contribute to the instanton interaction. The 
interaction between instantons is expected to suppress overlapping instantons and to result 
in a dynamical self-consistent cut-off. 



Some consequences of this picture have been discussed in [[□], [T2|], based on certain 
assumptions about the repulsive instanton interactions. In [|1^] a temporary infrared cut- 
off was introduced by means of a finite space-time volume V. The large V limit was 
then considered with the help of renormalization group arguments. This led to some 
general results independent of the specific form of the repulsive instanton interactions. In 
particular, a finite renormalization factor 

b = liiV (6) 
appears in some quantities, e.g. in the trace anomaly jll|| , correcting inconsistencies with 
the renormalization group. The same factor is conjectured to multiply the instanton 
singularities in the Borel plane, which then coincide with the infrared renormalons. The 
conclusions were supported by considering a model of the instanton ensemble, where the 
repulsive interaction is approximated by a hard core. 

The theory of instanton ensembles with a dynamical size cut-off has been developed 
further by Shuryak [|13[ 
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1C] for a review. In his model of an "instanton liquid" various 



observables have been calculated and related to hadronic phenomenology. 

In connection with the dynamical cut-off the distribution of instanton sizes is of central 
importance. The size-distribution reflects the way in which large instantons are suppressed 
and thus gives information about the instanton interactions. In recent years it has been 



studied by means of lattice Monte Carlo calculations by different groups [14], [15|, [Tj| [17 
For small sizes the distribution is predicted to be 



n{p) ~ p 
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(7) 



by the dilute gas approximation as well as by the "instanton liquid model" , in accordance 
with Eq. (|5]). For large sizes p, where the dynamical cut-off is in effect, not much is known 
about the distribution. There are arguments O, [L8|, |19[ in favour of a suppression like 



n(p) ~ exp(— cp p 



with 



P 



2. 



In this article we investigate the distribution of instanton sizes in a model |TJ] where 
the instanton interactions are approximated by a repulsive hard core of variable size. 



Although this approximation appears to be crude, the general features of the instanton 
ensemble with a dynamical cut-off are present. In particular, using analytical and numer- 
ical methods we calculate the asymptotic behaviour for small and for large sizes p and 
compare them with results from Monte Carlo simulations of lattice gauge theory. More 
details can be found in |30| . 

In 1 it has been conjectured that the distribution n(p) is affected by the finite 
renormalization factor 4/6 in such a way that for small p asymptotically 



n(p) 
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(9) 



Using the simplified model, we shall show below that this conjecture is wrong and that 
instead the semiclassical result (0) holds. 



1.2 Simplified model for ensembles of instantons 



Consider an ensemble of instantons in d space-time dimensions. In the spirit of ]TI| we 
introduce a finite volume V and study the approach to the thermodynamic limit V — > oo. 
In the sector with instanton number K the partition function is written as 

iK „ K K 



Z K (V) = ^- [Udz^Up^e 

n - J i=l j=l 



b-d-l-U({ ak },{ Pk }) 



(10) 



where the instanton positions and radii are denoted {aj,pj}. C is a constant, whose 
numerical value is unimportant here, b = UN/ 3 for SU(iV) Yang-Mills theory, and U 
represents the interaction potential between instantons. Distances are measured in units 
of A" 1 . 

In our simplified model the repulsive potential is approximated by a hard core poten- 
tial. The radius of an instanton core varies proportional to the size p of the instanton. In 
a finite volume V this means 
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is the volume of the unit sphere in d dimensions. The parameter r specifies the effective 
volume rp d : of an instanton and is of the order of V\. 



We introduce a reduced distribution by 

qK n K K-l K 



such that 



Z r K ed (V, PK) = -^fU da, I] d Pj II P b f d ^(M, {Pi}) , (13) 
n - J t=i i=i fc=i 

MV) = JdpZ™ d (V,p). (14) 

For the total system with variable instanton number the grand canonical partition 
function is 

oo 

Z(V) = Y, z k{V), where Z Q (V) = 1 . (15) 

K=0 

We do not distinguish between instantons and anti-instantons in this model. In this way 
we neglect aspects of the interactions which differ between instantons and anti-instantons, 
but we do not expect that they play a significant role for our considerations. 
The probability distribution of instanton numbers is given by 

In order to define the probability distribution of instanton sizes one has to specify 
how the sizes are sampled. The definition should be made in such a way that it is 
compatible with the Monte Carlo calculations to be discussed later. In the Monte Carlo 
runs configurations with a variable number of instantons are produced. A configuration 
of K instantons contributes K entries to the total histogram of instanton sizes. Therefore 
it has a relative weight proportional to K. Correspondingly the probability distribution 
in the iT-instanton sector is normalized to K [^TJ : 

«*<™-*w- (17) 

In the total ensemble the sizes are then distributed according to 

oo 

n(V,p) = J2 v K(V)n K (V,p) (18) 

oo 

YKZ r K ed (V,p) (19) 



K=l 

V KZ reil < 

Z(V) 



K=\ 

with 

PCX) 

/ n(V,p)dp=(K) v . (20) 
Jo 

As the expectation value of the instanton number grows linearly with the volume V one 
is interested in the rescaled distribution 

n(V, P ) = ^. (21) 

In the following sections we study the properties of n(V,p), and its thermodynamic limit 
n(p), respectively, utilizing analytical approaches as well as Monte Carlo- methods. 
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2 The one-dimensional instanton gas 



In d = 1 dimensions the model can be solved exactly in the thermodynamic limit. This 
case illustrates some general features and can serve as a testing ground for approximations 
used in higher dimensions. Therefore we shall discuss these results before we turn to the 
consideration of other dimensions. 

In the one-dimensional case the canonical partition function for a system of spatial 
length L can be written as 

(jK , K K 



Zk(L) = —JI[ dadPi II P b f 2 ®(M, { P k}) (22) 



i=l j=l 



with 



1.) -a,j\ > (f J (pi + pj) and 
©(R},{p*}) =1, if { 2.) |oi|<±L Vi, and 



3.) 0< Pi <± Vz 
©({^i}) {Pi}) = , else. 

This represents a system of rods with variable lengths on a line. A pictorial representation 
is given in Fig. |l[ 




rpx rp 2 rp 3 rp 4 

Figure 1: One-dimensional instanton gas, K = 4 



Integration over the instanton positions {a^} yields the effective free volume of K 
indistinguishable particles on the line L: 

K 



J f[dai 6(H), { Pfc }) = &({Pk}) (l-t£ P ^ 



with 



j=l \ 1=1 
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(23) 



e({p fe })H J' U^ iPj - L ^ 



as can be shown by induction. With this result the partition function reads 

r K r- r--T K P K ( K \ K 



The integrations over the pj can be carried out successively employing 

V(l-,)>d,= r (" + 1 ^ t V>, (26) 

and one obtains for the size distribution 

£g? =1 ATC^-Hr^ - 1) + 2))-y- 2 (L - rpff^+i 

nl ,p) E? =0 7 K (r(^ + 1))- 1 ^ 1 1 g 

with 

CT(6-1) , . 

7= ^ZI J - (28) 

The asymptotic behaviour for small instanton sizes is given by the power law 

n(L,p) ~ const, p 6 ~ 2 for p^O (29) 

in the sense of 

b-2 = hm H ^\ p)) . (30) 
p->o ln(p) 

In order to discuss the behaviour for large p we take the infinite volume limit 

n(L, p) L ^ n(p) . (31) 

The grand canonical sums can be evaluated by replacing them by integrals over K and 
performing a saddle point approximation, which becomes exact in the large- L limit. The 
result is 

n(p) = j p"- 2 *- cp , (32) 

with 

c= (CT{b-2)rf* . (33) 

In addition to the power law with exponent b — 2 = b — d — 1 we recognize an exponential 
suppression of large instanton sizes. The exponent p in the exponential, cp. Eq. (|8]), is 
equal to 1 in one dimension. Our next aim is to see how these results generalize to higher 
dimensions d. 



3 The general instanton gas 

In higher dimensions, d > 1, it is not possible to derive closed expressions for the partition 
functions or the size distributions. The main difficulty is that the integrations over the 
instanton positions and the radii cannot be decoupled. In particular we have to use 
approximations for the effective free volume of sets of instantons in higher dimensions. 
Nevertheless one can obtain approximate expressions for n(V, p) and derive its asymptotic 
behaviour for small p. 
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The canonical i^-instanton partition function Zk(V) is written as 

qK ~v n K K 

1,v / IT d&idp. 

r- 



Z K {V) = ^ f V dv f n dmdpt fl P] 6({a fc }, {p k }) S(v -rJ2pf), (34) 
■ J »=l 1=1 i 



where 



and 



a = 6 — ci — 1 , (35) 



rY.Pt <V (36) 



v 

i 

is the total effective volume of instantons. The maximal accessible volume V < V accounts 
for the fact that spheres in dimensions d > 1 cannot completely occupy a given volume 
V. 

An approximative decoupling of the integrations can be achieved through the obser- 
vation that for given v — r J2j Pj the product 

K 

IK 

3=1 

develops a sharp maximum at 

Pj = (jfr) =P°> j£{h-,K} (37) 

in the thermodynamic limit. The main idea is then to perform a saddle point approxima- 
tion for the p-integrations near this sharp maximum. The integrations over the positions 
aj then correspond to a gas of hard spheres with equal radii po- With 

Pj = Po + $j (38) 

and 

6(v -rfrf) » Si-rdpt 1 = ^ f<k e*^'' (39) 

the p-integrations can be solved straightforwardly. 

In a similar way the reduced partition function Z^ d (V, p) can be evaluated. For a given 
radius p of the K th instanton one assigns an effective volume V — rp d to the remaining 
K — 1 instantons and performs the saddle point approximation for the integrations over 
p±, . . . , pk-i in terms of Gaussian integrals. 

The integral over the instanton positions 

J J] da* e({ aj }, {p }) (40) 

3=1 

can be estimated with the help of geometrical considerations p2| , pi , 23], [1^, 0. We use 
an approximation of the form 

W) « ^(V - ^e//) K , (41) 

3=1 



where 

Ve f f = h(d)v, (42) 

and h(d) = V/V measures the inverse filling fraction of spheres in a volume V. In the 
one-dimensional case we get h(l) = 1 because a given length can be completely filled with 
rods. For d > 1 we consider h(d) as a parameter. Lower and upper bounds are given by 
1 < h(d) < 2 d - x . For more details on this point see gl], ||, |0 . 

Using these expressions we get by some lengthy but straightforward calculations for 
the canonical partition functions 

. is 

l)K) V? K , (43) 




T((3K + 1 



and for the reduced ones 
Z r K ed (V,p)tt 




K-l 



2tt d \ V a ((K - l)Th(d)y~ l ) F(J3(K-l) + 2 
where 

V = V- h(d)rp d , P=- r - (45) 

d 

The next step is to perform the grand-canonical sums over the instanton numbers 
K. As in the one-dimensional case, the sums can be evaluated in the large volume 
limit by replacing them by integrals which are calculated by means of the saddle point 
method. The error of this approximation vanishes in the thermodynamic limit. For the 
size distribution we get in this way 

(Id 

n(p) = P^" 1 exp(-cp rf ) , (46) 

with 




c *= c \l b T A 1 ,- 1 ) e~^h(d)r. (47) 



For b ^> d this takes the form 

C 3 =Cr(^-l)tH h(d)r, (48) 

which agrees with the one- dimensional result. 

The expression for n(p) is consistent with the general expectation mentioned in the 
introduction: for small p it grows powerlike with an exponent a — b — d— 1, and for large 
p this power-law is combined with an exponential decrease. 

Although the canonical partition functions are dominated by configurations where the 
instantons are densely packed, the exponent a agrees with the one of the semiclassical 
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dilute gas approximation. This result does not depend on the details of our approxima- 
tions and follows from the general structure of the occuring terms in the grand-canonical 



sums. The conjecture, made in |12| , that due to the denseness of instantons the small-p 
behaviour of n(p) gets modified, is therefore wrong. 

On the other hand, the value of the exponent p = d in the exponential decay at 
large p should be considered with reservations, because it depends on the saddle point 
approximations which have been made. In d — 1 dimensions it is correct, but we would 
not be surprised, if in higher dimensions the true value would differ from d. In order to 
get more insight into this question and to get an idea of the quality of the approximations 
being made so far, we have also studied the instanton gas by grand canonical Monte Carlo 
simulations. 



4 Monte Carlo simulations 

Usually Monte Carlo calculations are done in the canonical ensemble. In our case the 
particle number has to change and it is necessary to simulate a grand canonical ensemble. 
Simulations of grand canonical systems are not very common. They are rarely discussed 
in the literature and some important details remain unclear. Therefore it appears appro- 
priate to describe the algorithm we have used in our calculations. For related work on 



this topic we refer to [24, 25, 26, 27, 28, g| 



4.1 Grand canonical Monte Carlo algorithms 

A stochastic process, which is realized in a Monte Carlo simulation, is specified by a 
transition matrix W(X, Y), where X and Y denote states of the system. For the purpose 
of a simulation W is usually decomposed as a product of two factors: u(X, Y) represents 
a proposal probability for a transition from X to Y, and axy denotes the corresponding 
acceptance probability. In addition to normalization and ergodicity one has to require 
stationarity, which is often fulfilled by demanding the stronger Metropolis condition of 
detailed balance: 

' uj(Y,X)P(Y)\ 

Here P is the probability distribution, which we want to generate as the stationary dis- 
tribution of the underlying stochastic process. In our context it is given by 



P K (V; a 1 ,...,* K ,Pi,...,PK) = — ^ ft PJ (K'}> iPi}) ■ ( 5 °) 

The states X and Y are characterized by the instanton number K combined with the set 
of coordinates {a.j,pj}. 

In canonical algorithms u(X, Y) is usually chosen to be symmetric so that it is omitted 
in ( f49"D without further comments. This is not possible in a grand canonical ensemble, 
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where one has to consider transitions that change the instanton number. Independent 
of the choice of uj(X,Y) there will be additional volume factors in a X y for processes 
that do not conserve the instanton number. This results from the asymmetry in particle 
creation and destruction. If an instanton is created one has to specify a probability for 
the generation of its new coordinates, On the other hand, in the process of removing an 
instanton such a probability does not occur. In our case, for the space-coordinates as well 
as for the radii we choose a uniform distribution within the allowed volume. 

In the algorithm three different kinds of steps occur with equal probability: creation, 
destruction and movement of an instanton. With the shortcut notation 



1. 
2. 

3. 



> 



(Pi + Pj) Vi,j, and 



< \V d V/x, i, and 



we have chosen the following transition rules, where x denotes a random number between 
and 1. 



Creation: 

The creation of a new instanton with number K + 1 and coordinates (a', p') is 
proposed, and 

Y, CV ^ K v l{f p' a > x and * 

X, CV 2(K+l) ' d P^ < X 01 n0t * • 

• Destruction: 

The destruction of an instanton with randomly chosen number j is proposed, and 



X 



X 



Y, 
X, 



2K 



CV(Vv!/T)d 



r P j a >x 



IK -a ^ 

X Pi <~ X . 



CV(VV!/T)d 



• Movement: 

A movement of a randomly chosen instanton in a volume element [— 5 a , 5 a ] d x [S p , S p ] 
around the original coordinates is proposed, and 



X 



Y, i-^j >x and * 

X, i^-j < x or not * . 



The simulations were started with the empty configuration (K = 0). Measuring was 
started after the instanton number reached saturation. 
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4.2 Simulation results 

In the case of d = 1 dimensions the available exact result (^) provides a useful check 
of the Monte Carlo calculations. In Fig. |2| Monte Carlo data for n(p) in d = 1 are 
compared with the exact formula. The size L has been chosen large enough such that 
finite L effects are negligible. Obviously the Monte Carlo data agree very well with the 
theoretical predictions, and the thermodynamic limit has been approached sufficiently 

0.12 

0.1 

0.08 

n(p) 0.06 

0.04 

0.02 



012345678 

P 

Figure 2: The size distribution n(p) in d = 1 from a Monte Carlo simulation with C = 1, 
t = 2, a = 1, V = L = 2000 in comparison with the predictions of formula (|32|). 

With this check on the Monte Carlo algorithm we proceed to the more interesting case 
of four space-time dimensions (d = 4). In order to compare the Monte Carlo data with 
the outcome of our analytical approximations, Eq. fl4"E|), we have to make assumptions 
concerning the parameter h(d) = h{4) that describes the ability of instantons to fill a given 
volume. We consider three choices, namely the lower bound hi = 1, the upper bound 
h 2 = 2 4_1 = 8 and their geometric mean h g = 2y/2. The parameter r is taken to be 
equal to v±. For the volume we chose V = 15 4 . This is based on simulations in different 
volumes, which showed that in this case the thermodynamic limit was approximately 
reached within the errors of the simulation. The parameter a is taken to be a = 7/3, 
which is the value for SU(2) gauge theory in 4 dimensions. 

Fig. shows the Monte Carlo data in comparison with the analytical approximation. 
Near the maximum of the distribution the approximation qualitatively reproduces the 
Monte Carlo results. Furthermore, the growth of the distribution for small instanton 
radii according to a power law with exponent a can be confirmed, as is shown in Fig. [|. 
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Figure 3: The size distribution n(p) in d — 4 dimensions from a Monte Carlo simulation 
with C = 1, a = 7/3 (SU(2)), V = 15 4 in comparison with the predictions of formula 
(H) for h(A) = 1, 2^2, and 8. 
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Figure 4: The size distribution n(p) in d — 4 dimensions from a Monte Carlo simulation 
with C = 1, a = 7/3 (SU(2)), V = 15 4 in comparison with the predictions of formula 
for h(4) = 1, 2^2, and 8, plotted on a double logarithmic scale. 
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In order to study the behaviour of n(p) for large p we considered the ratio 

n(p) 



F(P) 



Inspired by the theoretical results we tried fits of the form 

Ffit(p) = aexp(-cpP) . 



(51) 



(52) 



The parameter a was obtained by extrapolating F(p) to small p. The fit with parameters 
c and p was then obtained using the Marquardt-Levenberg-algorithm. We performed fits 
for various choices of the model parameters C and a. In agreement with the theoretical 
results they showed that c depends on a, while p is nearly independent of it. 

The main interest is in the exponent p. We present the results for the parameter set 
C — 1, a — 7/3, V = 15 4 , because this value of a is relevant for gauge theory with gauge 
group SU(2). For a = 6, the SU(3) case, the results for p are the same within the present 
errors. We find a « 0.89, and the fit leads to c = 3.3 ± 0.2 and p — 1.9 ± 0.2. In Fig. | 
the result of a fit in the interval [0, 2.25] is shown. 



F(p) 




Figure 5: F(p) = n(p)/p a in d — 4 dimensions from a Monte Carlo simulation with 
C = 1, a = 7/3 (SU(2)), V = 15 4 in comparison with the fit F fit (p) with a = 0.89, 
c = 3.24, p = 1.92, plotted on a logarithmic scale. 



This Gaussian (p = 2) decay of the probability distribution has already been predicted 
by some authors under various assumptions JTTJ, [18] . A recent approach based on an idea 
of dual superconductivity |19|] also leads to the prediction p = 2. Furthermore good 
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agreement with the 577(3) lattice gauge theory calculations of Hasenfratz et al. was 
found. 

The exponent p = 2 differs from the one predicted by our approximate analytical 
calculation, p = d. The saddle point approximation being made is, however, not beyond 
any doubt. In that case the exponent originates from the effective excluded volume being 
proportional to p d at the considered saddle point. This would also coincide with the 
intuitive expectation based on the following picture. In the presence of a large instanton 
of size p the remaining ones are excluded from a volume ~ p d . If they behave like a dilute 
gas, one would expect that the excluded volume yields a suppression factor oc exp(— cp d ). 
The instanton ensemble in the effective remaining volume is, however, dominated by 
dense configurations, as the analytical calculation shows. Therefore the intuitive picture 



should be considered with reservations. Indeed, the calculations of |TlJ take into account 
excluded volume effects in the framework of the theory of grand canonical pair distribution 
functions and, also employing certain approximations, arrive at p = 2. 

In recent years much effort has been devoted to lattice Monte Carlo calculations of 
properties of the instanton ensemble. There are still ambiguities due to smoothing proce- 
dures and only data with little statistics are yet available. Nevertheless some quantitative 
statements have been given. Concerning the size distribution for small p, lattice calcu- 
lations appear to support the power law © rather than (g). A nice plot, using data 
of [|I6|1 , can be found in |3T| . For the large-p distribution, de Forcrand et al. predict an 
exponential decrease with p = 3 ± 1 from their SU(2) lattice data jn|. In contrast to 
this, Smith and Teper conclude form their 577(3) simulations a decay according to p~^ 
with f « 10 . . . 12 |T§ . 

We have studied the distribution of instanton sizes p in the framework of a model, 
where instanton interactions are approximated by a hard core potential with variable 
radius. This model incorporates the basic features of a dynamical cut-off on large instan- 
ton sizes. In the one-dimensional case an exact formula can be derived, which yields a 
power-like growth ~ p a for small radii p and an exponential decay for large p. 

In four space-time dimensions we employed analytical approximations as well as Monte 
Carlo simulations. The theoretical calculations generalize the one-dimensional results and 
give a power-like behaviour for small p. For large radii p they overestimate the decay which 
is found in the Monte Carlo data. Fits to the numerical Monte Carlo results suggest a 
behaviour like 

n(p) exp(-cp 2 ) , (53) 

in agreement with some other work on gauge theories. 

The results indicate that our simplified model reproduces the main features of in- 
stanton ensembles with a dynamical infrared cut-off. Definite results about properties of 
instanton ensembles can of course only be expected from future Monte Carlo calculations 
of lattice gauge theories. 
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